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We present a theoretical framework for a quantitative understanding of the full doping dependence 
of the optical spectra of the cuprates. In accord with experimental observations, the computed 
spectra show how the high-energy Mott features continue to persist in the overdoped regime even 
after the mid-infrared (MIR) peak originating from the pseudogap has collapsed in a quantum 
critical point. In this way, we reconcile the opposing tendencies of the MIR and Mott peaks to shift 
in opposite directions in the optical spectra with increasing doping. The competition between the 
pseudogap and the Mott gap also results in rapid loss of spectral weight in the high energy region 
with doping. 

PACS numbers: 74.20.-z, 74.20.Mn,74.25.Gz, 74.25. Jb 



I. INTRODUCTION 

Most experiments on the cuprates show that when the 
Mott insulating state is doped with electrons or holes the 
gap feature in the spectrum moves to lower energies and 
collapses at a quantum critical point, consistent with the 
notion that electron correlation effects in the Mott state 
weaken systematically with doping In sharp contrast, 
optical experiments reveal a substantially more complex 
picture in that the doping of the Mott insulator induces a 
new mid-infrared (MIR) feature^—, which collapses with 
doping much as seen in other experiments, but at the 
same time a high-energy Mott feature persists even into 
the overdoped system with its spectral weight shifting 
to the MIR feature with increasing doping In fact, 
the Mott feature not only persists in the optical experi- 
ments, but it moves to higher energies with doping, sug- 
gesting that strong electron correlation effects continue 
to play an important role in the cuprates at all dopings. 
Modelling and understanding the optical spectra of the 
cuprates thus becomes of key importance in unravelling 
the routes by which the Mott insulator turns itself into a 
superconductor, a transition that remains poorly under- 
stood. 

The question of correlations has been traditionally 
framed in terms of a Slater picture of itinerant electrons 
where an insulator forms a gap via the development of 
a long range magnetic order, or a Mott picture in which 
the metal-insulator transition is driven by a local condi- 
tion of no double occupancy. Although the Slater or the 
Mott picture is often invoked, many materials lie in the 
crossover regime where the electronic states are not well 
described as being either fully itinerant or fully local- 
ized. Such materials are often endowed with unique and 
exotic properties and are of great current interest from 
the viewpoint of fundamental physics as well as potential 
for applications. The strength of electron correlations in 
the cuprates has been a matter of considerable debate 
since the discovery of these fascinating materials 

We have obtained the optical spectra within the frame- 
work of a one band Hubbard model where the self-energy 



is obtained self-consistently in an intermediate coupling 
scenario to account for spin and charge fluctuations via a 
computation of the susceptibility over the entire doping 
range. We emphasize that the present model essentially 
does not involve any free parametrization in that the bare 
LDA dispersion is taken to be the CuC>2 band (in the 
tight-binding form) and a doping-independent bare U is 
chosen to reproduce the experimental charge-transfer gap 
at half-filling. The screened U at various doping levels 
is then computed due to charge fluctuation o 15 ' 16 , which 
self-consistently induce a low-energy magnetic order in 
the in-gap states resulting in a pseudogap in the under- 
doped system. Our computations reproduce quantita- 
tively the experimentally observed doping evolution of 
not only the shifts in the positions of the MIR and Mott 
features, but also the rapid transfer of spectral weight 
from the high to the low energy region with increas- 
ing doping. The theoretical shifts of the MIR feature 
with doping are in accord with angle-resolved photoe- 
mission spectroscopy (ARPES)i£ and other spectroscopic 
probes. We emphasize that existing approaches invoking 
either strong or weak coupling scenarios fail to capture 
the essence of the doping evolution of the optical spectra. 



II. OPTICAL SPECTRA IN ELECTRON AND 
HOLE DOPED CUPRATES 

Fig. 1 shows one of our key results. The computed 
evolution of the optical conductivity <t(uj) is seen to be 
in excellent accord with measurements on electron-doped 
Nd 2 _ x Ce :[: Cu04 (NCCO) and hole doped La 2 - a; Sr a; Cu04 
(LSCO) 1:L i 13 . All the spectra show a nearly isosbetic or 
equal absorption point near 1.3 eV (1 eV for LSCO) 
[black vertical line], consistent with the experimental 
behavior™. The doping evolution is completely differ- 
ent on opposite sides of this isosbetic point. Above this 
point, the spectrum is dominated by a broad hump fea- 
ture associated with the Mott gap. At half-filling, only 
this feature is present and the calculated optical spec- 
tra show an insulating gap whose energy, structure, and 
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FIG. 1: (color online) Comparison between calculated optical 
spectra and experiments in NCCO in (a) and LSCO in (b). 
Experimental results for NCCO are taken from Ref.[l3l. except 
for x — 0.17, which including LSCO data, are from Ref. [Til . 
The x = 0.17 dataset for NCCO includes a background sub- 
traction to match the former dataset^. Dashed line gives the 
background contribution added to the theoretical spectrum 
at x = [Ref. [3|. 



intensity match remarkably with measurement a 1 1 1 13 . As 
doping increases, the high energy peak shifts to higher en- 
ergy and broadens and its spectral weight is transferred 
to the Drude and MIR peaks. The MIR peak shifts to 
lower energy with doping and gradually sharpens. Note 
that in both samples at the highest doping, when the 
MIR peak collapses into the Drude peak, Mott-gap fea- 
tures still persist in the spectrum. The doping evolu- 
tion also shows similar behavior in other cuprates ^ 13 ' 19 . 
Notably, our calculations also describe the anomalous 
<7 ~ 1/w-dependence found in most cuprates associated 
with magnetic scattering. 



as is also the case in variational cluster calculations 21 . 
The highest and lowest bands are an incoherent residue of 
the undressed bands, which we will refer to as upper and 
lower Hubbard bands. The two inner bands are coherent 
in-gap states split by a spin density wave induced AFM 
gap into upper and lower magnetic bands. The in-gap 
states and the high energy Hubbard bands are separated 
by high-energy kinks ('waterfalls') predominantly associ- 
ated with magnetic excitations, as observed universally 
in all cuprates by ARPES 22 ~— , and found in quantum 
monte carlo (QMG)^a and variational 21 calculations. 

The corresponding optical spectrum in Fig. 2(b) con- 
sists of three main regions marked by shading of dif- 
ferent color: (1) The low frequency Drude region for 
u < 30meV; (2) The MIR region; and, (3) the high en- 
ergy Mott gap region for ui > 1.5eV. We concentrate 
here on the MIR and Mott-gap regions, and return be- 
low to comment on the Drude region. The interband 
optical absorption is proportional to the joint density of 
states (JDOS), so that at each energy we can construct 
a 'source map' showing the filled states which make a 
strong contribution to the transition, and a 'sink map' 
showing the contribution of the corresponding empty 
states. Figs. 2(d) and (e) show the states responsible 
for the optical transitions along the high symmetry lines 
at representative photon energies of ui = 0.5 eV near MIR 
peak and uj — 2 eV around the high energy peak (vertical 
bars in Fig. 2(b)). At cj = 0.5 eV, the transitions are con- 
fined within the in-gap states only, whereas at u = 2 eV, 
the optical subbands involve predominantly the incoher- 
ent region. The depletion in the optical spectral weight 
near the isosbetic point in Fig. 1 and in Fig. 2(b) is thus 
associated with the 'waterfall' region marked by arrows 
in the ARPES spectrum in Fig. 2(a)±£. 

The total optical spectral weight obtained in Fig. 2(b) 
is the integral of the JDOS times the band velocity. The 
role of the latter factor is explicated in Fig. 2(c), where 
contributions to a are plotted as a function of photon en- 
ergy uj and momentum k. Although the source-sink map 
shows a symmetry about the (tt, 0) point, the quasipar- 
ticle velocity is low along the T — » (ir, 0) direction [cyan 
solid line in Fig. 2(c)], leading to greatly reduced spectral 
intensity associated with those regions. In contrast, the 
large quasiparticle velocity in the other two directions is 
responsible for two distinctive intense streaks (labelled d 
and e). The MIR peak is clearly dominated by the antin- 
odal quasiparticles, whereas the high energy hump stems 
from a wider fc-range. 



A. Origin of optical dichotomy 

Fig. 2 helps delineate the microscopic origin of these 
features by comparing the spectral intensities relevant for 
ARPES and optical spectra at a representative doping of 
x = 0.10 for NCCO; similar analysis for LSCO is not 
shown for brevity. In the computed ARPES spectrum in 
Fig. 2(a), the underlying LDA dispersion is clearly visi- 
ble, but the spectral weight has split into four subbands, 



B. Competition between Pseudogap amd Mott gap 

Fig. 3 summarizes our results and elaborates upon the 
competition between the collapse of the pseudogap in 
the coherent bands at a critical doping, and the per- 
sistence of the Mott gap in the incoherent bands. The 
gaps extracted from the optical spectra are now com- 
pared with other experimental data for various electron- 
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FIG. 2: (Color online) Connection between optical spectra and spectral intensity maps in NCCO. (a) Computed spectral 
intensity (on log scale) as a function of w along the high symmetry lines at a representative doping of x — 0.10. (b) Calculated 
optical spectra are compared with the corresponding experimental results^. Green solid line gives the spectrum calculated with 
coherent bands only, i.e. without the self-energy. Dashed green line is based on using LDA bands with self-energy correction 
but without the pseudogap. The shaded regions of different colors (left to right) approximately mark Drude, MIR and high 
energy regions. Two vertical black lines indicate the photon energies at which frames (d) and (e) are calculated, (d), (e) give 
source and sink maps corresponding to the optical transitions at fixed photon energy cj along the high symmetry directions. 
White vertical arrows indicate the photon energy connecting source and sink points involved in a particular transition. White 
shaded region in (e) highlights the incoherent part of the spectral weight, not visible in (d). (c) Optical spectrum as a function 
of photon energy is plotted along the high symmetry lines. Cyan line gives the band velocity as a function of k (arbitrary 
units). The two white arrows indicate the two photon energies at which (d) and (e) are calculated. 



and hole-doped cuprates. The pseudogap seen in the op- 
tical spectrum corresponds to the true AFM gap, and 
its doping evolution is in good agreement with opti- 
cal, ARPES^ i 26 ' 31 and resonant inelastic x-ray scattering 
(RIXS)22. results. Both ARPES and MIR data predict a 
QCP near x = 0.1 7 32 ' 33 consistent with the present the- 
oretical predictions. In contrast, the Mott gap shows the 
opposite doping dependence to the pseudogap, increas- 
ing slowly with doping. Although the Mott gap does not 
show a real QCP, it rapidly loses intensity with doping. 

Similar results are found in hole doped LSCOii, 
and other layered cuprates, including Bi 2 Sr 2 CaCu 2 08+,5 
(B12212) 9 , YBa 2 Cu 3 6+:!; (YBCO)£, Tl 2 Ba 2 Cu0 6+(5 
(TBCO), and HgBa 2 Ca2Cu 3 8+< 5 (HBCO)^ as well 
as in x-ray absorption spectroscopy^ and QMC 
computations^ Figure 3 also includes data on the MIR 
gap for LSCO, Bi2212 and YBCO and the pseudogap ob- 
tained from a wide range of experimental data, including 
a recent detailed survey^ (see Fig. 3 caption). In the 
intermediate and high doping regimes, there is very good 
agreement between the MIR and other measures of the 



pseudogap, while the MIR data provides important evi- 
dence for the rapid growth of the pseudogap in the deeply 
underdoped regime, consistent with theory. We note that 
while all the hole-doped cuprates seem to have a very 
similar doping dependence of the pseudogap, there are 
subtle differences with electron doped cuprates, includ- 
ing the steepness of the rise at low doping and the exact 
position of the QCP i 28 i 35 The former difference seems to 
be correlated with the doping dependence of the screened 
Hubbard J7's, also shown in Fig. 3. Since the doping de- 
pendence is mainly a screening effect ) 15 ' 16 the difference 
is probably due to the strong screening associated with 
the van-Hove singularity (VHS) on the hole-doped side. 

While a strong case can be made that in electron- 
doped cuprates the pseudogap is associated with a co- 
existing (jr, 7r) AFM order, the nature of the pseudogap 
in hole doped cuprates is not well understood. There is 
growing consensus that it originates from some form of 
density-wave like competing order (which could include 
coupling to phonons)^ 6 - - — , and may actually also involve 
several competing modes, again due to proximity to the 
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FIG. 3: (Color online) Contrasting doping dependences of 
the MIR and Mott gaps in the cuprates. Present theoretical 
results are compared to optical dat a 10 ' 11 ' 13 and pseudogap 
data measured by ARPES^ 6 - and RIX$2£ for electron doped 
NCCO and many other probes for hole doped cuprates. The 
red squares are the pseudogap for four different layered com- 
pounds [Bi2212, YBCO, TBCO and HBCO] as they are ob- 
tained from ARPES, tunneling, Raman, Andreev reflection, 
and heat capacity experimental data (reproduced from a de- 
tailed survey, Ref. [28T h The red diamond and star symbols 
represent independent MIR and green diamond and star give 
FIR gaps measured simultaneosly in optical spectra of Bi2212 
and YBCQ22, while the green circles represent a Coulomb gap 
seen m ARPES 3 ^ [energy values are obtained assuming a ratio 
2A/fcsT = 4]. In all cases, open symbols of same colors are 
the corresponding theory which are obtained self-consistently 
at each doping. Shown also are the computed screened value 
of U (solid lines) 15 ' 16 used in these calculations. The dashed 
lines are guides to the eye for the Mott gap (gold), the pseu- 
dogaps for electron (blue) and hole (red) doping, and the FIR 
and Coulomb gaps (green). 



VHS^. However, in our earlier study in Ref. [H, which 
compared several possible competing order phases, we 
demonstrated that the shape and doping dependence of 
the pseudogap is remarkably insensitive to the particular 
order so long as the competing order vanishes in a QCP 
near optimal doping. 

While our present calculations reproduce the gap 
structures of cuprates, particularly on energy scales > 
100 meV, additional effects can arise at lower energies, 
including coupling to phonons and impurities. In particu- 
lar, Lupi et alS^ find in several cuprates that in addition 
to the MIR feature there is an additional far-infrared 
(FIR) peak which appears near 10% doping and shifts 
to higher energies at lower doping (green diamond and 
star symbols in Fig. 3). This seems to be a disorder ef- 
fect, and has a similar doping dependence to the so-called 
Coulomb gap seen in ARPES (green circles) in the low- 
temperature region of Bi2212^ 

We emphasize that the dichotomy between the pseu- 
dogap and Mott gap and the presence of the QCP are 
robust features of cuprates as the computations do not 




FIG. 4: (Color online) Optical sum rule and spectral weight 
transfer, (a) Effective number of electrons N e ff(u>) calculated 
from the optical spectra in Fig. 1 is compared with experimen- 
tal data^ii 3 - on NCCO. Results for LSCO are similar and are 
not shown for brevity, (b) Low energy weight N e f j at oj — 30 
meV as an estimate of Drude weight, compared with experi- 
mental results on various cuprates^r-^. Green line gives the 
direct computation of the normalized Drude weight n/m*. 



involve any free parametrization except the bare doping- 
independent value of U = 1.7 cV at half-filling. We have 
computed the doping dependence of U due to the screen- 
ing effects of charge fluctuation o 15 ' 16 and obtained effec- 
tive U values are shown Fig. 3. 



III. OPTICAL SUM RULE 

Finally, we use the integrated optical spectral weight 
to illustrate the Mott gap collapse. The effective electron 
number (per Cu atom), N e ff(uj), can be defined in terms 
of the optical conductivity integrated up to an energy w. 



2m V 
ne 2 hN 



(1) 



where tuq and e are the free electron mass and charge, 
respectively, and N is the number of Cu-atoms in a cell 
of volume V. The results in Fig. 4(a) show how rapidly 
spectral weight shifts to low energies with doping, cor- 
rectly reproducing the experimental behavior as well as 
dynamical mean-field calculations — , but incompatible 
with strong coupling models (such as the t — J or U — > oo 
Hubbard model) where one assumes no double occupancy 
of Cu sites. The present intermediate coupling model of 
Mott gap collapse on the other hand properly captures 
these spectral weight transfers as a function of doping. 

Our model also predicts the Drude weight oc 
X)k rik / m k' wn i cn can he compared to an experimen- 
tal estimate, taken as N e ff at a characteristic energy 
lu = 30 meV in Fig. 4(b) for NCCOi2, and with results 
for Pr 2 _ a; Ce :r Cu04 (PCCO)i£ and LSCO±i. Even here 
the agreement is quite good. Interestingly, at low dop- 
ing the Drude weight increases almost linearly with x. 
Within our model, this simply reflects scaling with the 
area of the FS pockets in the pseudogap state. At op- 
timal doping, the weight shows a sharp jump (green ar- 
row) associated with the appearance of the hole pocket. 
This topological FS transition is an intrinsic feature of 
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the model and has been found in NCCO near x ~ 0.15 
by several experimental probes such as ARPES 2 , Hall 
effect, and superconducting penetration depth^ 2 . For 
LSCO (blue symbols), experiments^ show a similar lin- 
ear behavior of the Drude weight, which also corresponds 
to the doping dependence of the area of the FS pocket^. 
The peak around x ~ 0.2 corresponds to the doping of 
the Van Hove singularity. 

IV. INTERMEDIATE COUPLING MODEL 

We evaluate the self-energy S as a convolution over 
the green function G and the interaction W ~ U 2 \ 
(including the full spectrum of charge and spin fluctu- 
ations within the random phase approximation (RPA)) 

a3 17 i 42 d 3 ; 

£(k j( 7, lWn ) = i;y 2 Z^ry CT , CT , J ^ 

G(k + q, a', iu> n + w p )r(k, q, iu) n , w p )Im[xRp A (q, w p )]. 

(2) 

where a is the spin index and the prime over the q sum- 
mation means that the summation is restricted to the 
magnetic Brillouin zone. Here the spin degree of free- 
dom r\ a .a' takes the value of 2 for the transverse di- 
rection and 1 for both longitudinal and charge modes. 
In the underdoped region, the pseudogap is modelled 
by an antiferromagnetic (AFM) order parameter, result- 
ing in G, x an d E becoming 2x2 tensors^. We de- 
fine a total self-energy as E* = USf\ + E, where f\ 
is the Pauli matrix along the a;— direction and US is 
the AFM gap defined below. The self-energy E* con- 
tains essentially two energy scales: (i) it gives rise to 
the SDW with an additional renormalization of the over- 
all quasiparticle dispersions in the low energy region, 
and (ii) at higher energies it produces the Hubbard 
bands. We use a modified self-consistent scheme, referred 
to as quasiparticle— GW (QP— GVF)-scheme in which G 
and W are calculated from an approximate self-energy 
Sq(oj) = USf\ + (l — Z^ 1 ) oil, where the renormaliza- 
tion factor Z is adjusted self-consistently to match the 
self-energy E* at low energy 16 i 24 i 43 The vertex correction 
r(k, q, u),u)p) in Eq. [5] is taken as its first order approx- 
imation (Ward's identity) as r(k, q, uj,uj p ) — 1/Z. Since 
the ^-dependence of E is weak2i, we further simplify the 
calculation by assuming a fc-independent E, which we 
calculate at a representative point k = (7r/2,7r/2). 

Our model is nearly parameter free. We take the dis- 
persions directly from the LDA calculations, accurately 
fitted by a one band tight-binding model^, without any 
adjustment of the resulting parameters. The suscepti- 
bilities are calculated numerically, with no further ap- 
proximations, while the renormalization constant Z is 
determined self-consistently. We have performed realis- 
tic Kanamori screening calculations due to charge fluc- 
tuations to obtain a self-consistent value of U and order 



parameter <fi at each doping starting from a doping inde- 
pendent bare U = 1.7 eV which reproduces the charge 
transfer gap at half-fillingi 15 i 16 

The optical conductivity, cr(w), is calculated using 
the standard linear response theory using the above 
E— dressed single particle states. At high energies, we 
subtract off a linear-in-cj background from the experi- 
mental spectra associated with interband transitions to 
higher-lying bands not included in the present one-band 
calculations, consistent with supplementary Fig. 4 of 
Ref. [l3- While keeping fixed the calculated weights of 
the Drude peaks, we have phenomenologically broadened 
the spectra with additional constant impurity scattering 
rates, similar to the universal form found in experiments 
on NCCG±2, LSCOii, and Bi2212^. 



V. DISCUSSION AND CONCLUSION 

While the (tt, tt) AFM order is known to be a robust 
feature in electron doped cuprates, in the hole doped case 
the exact nature of the pseudogap is unknown, and is 
likely to be incommensurate and possibly associated with 
charge order. This is consistent with our earlier study of 
ARPES^, where we compared several possible compet- 
ing ordered phases and demonstrated that the shape and 
doping dependence of the pseudogap is remarkably in- 
sensitive to the particular order so long as the competing 
order vanishes in a quantum critical point near optimal 
doping. This conclusion should hold even more strongly 
for the optical spectrum. Note that we have included all 
the charge and spin fluctuations in our calculations - the 
particular leading instability will only lead to a small re- 
arrangement of spectral weight in the low energy regime, 
leaving most of the optical spectrum virtually unchanged. 

We note that recent Gutzwiller approximation (GA) + 
RPA calculations find a number of nesting instabilities in 
hole-doped cuprates, resulting in a number of competing, 
incommensurate density-wave orders^ Since the insta- 
bilities depend on the bare susceptibilities, they are lo- 
cated at virtually the same q- values in both the spin and 
charge (electron-phonon) channels. Thus, the resulting 
optical spectra should be very similar in both cases. 

Finally, we note that the same model can be used to 
explore the role of magnetic fluctuations as the pairing 
bosons for superconductivity^, While the study of Ref. 
l46l revealed a significant contribution from magnetic fluc- 
tuations, it was found that fluctuations in different fre- 
quency ranges could act cooperatively to greatly enhance 
T c . This suggests that phonons could play a similar co- 
operative role which could explain the isotope effect^ 

In summary, we have shown that our model frame- 
work explains the salient features of the optical spec- 
tra of cuprates and their evolution with doping at a 
quantitative level. The rapid loss of high energy spec- 
tral weight and the associated shift of the MIR peak to 
low energies in the optical spectra with increasing dop- 
ing reflects collapse of the pseudogap order [here taken 
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as a (w, 7r) AFM], and the presence of a QCP near op- 
timal doping in the coherent in-gap states. By con- 
trast, the Mott gap in the incoherent states persists 
at all dopings including, in particular, the overdoped 
regime. The aforementioned coherent and incoherent 
states are connected via a high-energy kink driven pre- 
dominately by magnetic excitations. In the magnetic ex- 
citation spectrum, the dominant excitations lie in the 
waterfall region, while the remnant of the low-energy 
magnetic resonance mode^S is less significant 24 . Our 
model self-energy scheme would provide a tangible basis 
for modelling other spectroscopies (e.g., angle- resolved 
photoemissio n 17 i 43 , scanning-tunnelling^, inelastic light 
scattering 27 , positron-annihilation 47 ) of the cuprates and 



other materials. 



Acknowledgments 

This work was supported by the US Department 
of Energy, Basic Energy Sciences contract DE-FG02- 
07ER46352, and benefited from the allocation of su- 
percomputer time at NERSC and Northeastern Univer- 
sity's Advanced Scientific Computation Center (ASCC). 
RSM's work has been partially funded by the Marie Curie 
Grant PIIF-GA-2008-220790 SOQCS. 



1 A. Damascelfi, Zahid Hussain, and Zhi-Xun Shen, Rev. 
Mod. Phys. 75, 473 (2003). 

2 N.P. Armitage, F. Ronning, D. H. Lu, C. Kim, A. Damas- 
celli, K. M. Shen, D. L. Feng, H. Eisaki, Z.-X. Shen, P. K. 
Mang, N. Kaneko, M. Greven, Y. Onose, Y. Taguchi, and 
Y. Tokura, Phys. Rev. Lett. 88, 257001 (2002). 

3 M. Le Tacon, A. Sacuto, A. Georges, G. Kotliar, Y. Gallais, 
D. Colson, and A. Forget, Nature Physics 2, 537 (2006). 

4 K. Haule, and G. Kotliar, Europhys. Lett. 77, 27007 
(2007). 

5 S. Chakraborty, Dimitrios Galanakis, and Philip Phillips, 
Phys. Rev. B 78, 212504 (2008). 

6 Mishchenko, N. Nagaosa, Z.-X. Shen, G. De Filippis, V. 
Cataudella, T. P. Devereaux, C. Bernhard, K. W. Kim, 
and J. Zaanen, Phys. Rev. Lett. 100, 166401 (2008). 

7 L. Vidmar, J. Bonca, and S. Maekawa, Phys. Rev. B 79, 
125120 (2009). 

8 S. L. Cooper D. Reznik, A. Kotz, M. A. Karlow, R. Liu, M. 
V. Klein, W. C. Lee, J. Giapintzakis, D. M. Ginsberg, B. 
W. Veal and A. P. Paulikas, Phys. Rev. B 47, 8233 (1993). 

9 J. Hwang, T Timusk, and G D Gu, J. Phys.: Condens. 
Matter, 19, 125208 (2007). 

10 T. Arima, Y. Tokura, and S. Uchida, Phys. Rev. B 48, 
6597 (1993). 

11 S. Uchida, T. Ido, H. Takagi, T. Arima, Y. Tokura, and S. 
Tajima, Phys. Rev. B 43, 7942 (1991). 

12 Y. Onose, Y. Taguchi, K. Ishizaka, and Y. Tokura, Phys. 
Rev. Lett. 87, 217001 (2001). 

13 Y. Onose, Y. Taguchi, K. Ishizaka, and Y. Tokura, Phys. 
Rev. B 69, 024504 (2004). 

14 A. Comanac, Luca de Medici, Massimo Capone, and A. J. 
Millis, Nature Physics 4, 287 (2008). 

15 J. Kanamori, Prog. Theor. Phys. 30, 275 (1963). 

16 R.S. Markiewicz, and A. Bansil, Phys. Rev. B 75, 020508 
(2007). 

17 Susmita Basak, Tanmoy Das, Hsin Lin, J. Nieminen, M. 
Lindroos, R. S. Markiewicz, and A. Bansil, Phys. Rev. B 
80, 214520 (2009). 

18 Interestingly, the isosbetic point in the optical spectra as- 
sociated with the 'waterfall' features in ARPES occurs at 
lower energy in LSCO than NCCO, as does the high energy 
kink in ARPES spectral. 

19 McGuire, M. Windt, T. Startseva, T. Timusk, D. Colson, 
and V. Viallet-Guillen, Phys. Rev. B 62, 8711 (2000). 



20 C.T. Chen, F. Sette, Y. Ma, M. S. Hybertsen, E. B. 
Stechel, W. M. C. Foulkes, M. Schulter, S-W. Cheong, A. 
S. Cooper, L. W. Rupp, Jr., B. Batlogg, Y. L. Soo, Z. H. 
Ming, A. Krol, and Y. H. Kao, Phys. Rev. Lett. 66, 104 
(1991). 

21 C. Grober, R. Eder, and W. Hanke, Phys. Rev. B 62, 4336 
(2000). 

22 J. Graf, G.-H. Gweon, K. McElroy, S. Y. Zhou, C. Jozwiak, 
E. Rotenberg, A. Bill, T. Sasagawa, H. Eisaki, S. Uchida, 
H. Takagi, D.-H. Lee, and A. Lanzara, Phys. Rev. Lett. 
98, 067004 (2007). 

23 B. Moritz, F. Schmitt, W. Meevasana, S. Johnston, E. M. 
Motoyama, M. Greven, D. H. Lu, C. Kim, R. T. Scalettar, 
Z.-X. Shen, T. P. Devereaux, New Journal of Physics 11, 
093020 (2009). 

24 R.S. Markiewicz, S. Sahrakorpi, and A. Bansil, Phys. Rev. 
B 76 , 174514 (2007). 

25 A. Macridin, M. Jarrell, Thomas Maier, and D. J. 
Scalapino, Phys. Rev. Lett. 99, 237001 (2007). 

26 H. Matsui, T. Takahashi, T. Sato, K. Terashima, H. Ding, 
T. Uefuji, and K. Yamada, Phys. Rev. B 75, 224514 (2007). 

27 Y. W. Li, D. Qian, L. Wray, D. Hsieh, Y. Xia, Y. Kaga, 
T. Sasagawa, H. Takagi, R. S. Markiewicz, A. Bansil, H. 
Eisaki, S. Uchida, and M. Z. Hasan, Phys. Rev. B 78, 
073104 (2008). 

28 S. Hiiefner, M.A. Hossain, A. Damascelli, and GA. 
Sawatzky, Rep. Prog. Phys. 71, 062501 (2008). 

29 S. Lupi, D. Nicoletti, O. Limaj, L. Baldassarre, M. Or- 
tolani, P. Calvani, S. Ono, and Yoichi Ando, Preprint at 
larXiv:0905.0568l 

30 Z. -H. Pan, P. Richard, Y. -M. Xu, M. Neupane, P. Bishay, 
A. V. Fedorov, H. -Q. Luo, L. Fang, H. -H. Wen, Z. Wang, 
and H. Ding, Phys. Rev. B 79, 092507 (2009). 

31 C. Kusko, R. S. Markiewicz, M. Lindroos, and A. Bansil, 
Phys. Rev. B 66, 140513(R) (2002). 

32 Tanmoy Das, R.S. Markiewicz, and A. Bansil, Phys. Rev. 
Lett. 98, 197004 (2007). 

33 Tanmoy Das, R.S. Markiewicz, and A. Bansil, Phys. Rev. 
B 77, 134516 (2008). 

34 M. Jarrell, Th. Maier, M. H. Hettler3and A. N. Tahvil- 
darzadeh, Europhys. Lett. 56, 563 (2001). 

35 J. L. Tallon, and J. G. Storey, arXiv: 0908. 44301 

36 Jinho Lee, K. Fujita, K. McElroy, J. A. Slezak, M. Wang, 
Y. Aiura, H. Bando, M. Ishikado, T. Masui, J.-X. Zhu, A. 



7 



V. Balatsky, H. Eisaki, S. Uchida and J. C. Davis, Nature 
442, 546 (2006). 

A. Bussmann-Holder, H. Keller, A. R. Bishop, A. Simon, 
K. A. Mueller, larXiv:0809.288i1 

A. Macridin, and M. Jarrell, Phys. Rev. B 79, 104517 

(2009) . 

R.S. Markiewicz, J. Lorenzana, G. Seibold, and A. Ban- 
sil, Phys. Rev. B 81, 014509 (2010); R.S. Markiewicz, 
J. Lorenzana, and G. Seibold, Phys. Rev. B 81, 014510 

(2010) . 

Tanmoy Das, R.S. Markiewicz, and A. Bansil, Phys. Rev. 
B 74, 020506 (2006); J. Phys. Chem. Solids 69 2963 (2008). 
T. Yoshida, X J Zhou, D H Lu, Seiki Komiya, Yoichi Ando, 
H Eisaki, T Kakeshita, S Uchida, Z Hussain, Z-X Shen and 
A Fujimori, J. Phys.: Cond. Matt. 19, 125209 (2007). 



G. Vignale and M. R. Hedayati, Phys. Rev. B., 42, 786 
(1990). 

R. S. Markiewicz, Tanmoy Das, Susmita Basak, A. Ban- 
sil, International Workshop on Strong Correlations and 
Angle-Resolved Photoemission Spectroscopy, 2009, Zurich, 
Switzerland , to be published. 

J. R. Schrieffer, X. G. Wen, and S. C. Zhang, Phys. Rev. 
B. 39, 11663 (1989). 

R.S. Markiewicz, S. Sahrakorpi, M. Lindroos, Hsin Lin, 

and A. Bansil, Phys. Rev. B 72, 054519 (2005). 

R. S. Markiewicz, and A. Bansil, Phys. Rev. B 78, 134513 

(2008). 

L. A. Bansil, U. Welp, Y. Fang and K. G. Bailey, J. Phys. 
Chem. Solids 52, 1541 (1991). 



